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Abstract: We extract exact charged black- hole solutions with flat transverse sections in 
the framework of D-dimensional Maxwell-/ (T) gravity, and we analyze the singularities 
and horizons based on both torsion and curvature invariants. Interestingly enough, we find 
that in some particular solution subclasses there appear more singularities in the curvature 
scalars than in the torsion ones. This difference disappears in the uncharged case, or in the 
case where f(T) gravity becomes the usual linear-in-T teleparallel gravity, that is General 
Relativity. Curvature and torsion invariants behave very differently when matter fields are 
present, and thus f(R) gravity and f(T) gravity exhibit different features and cannot be 
directly re-casted each other. 
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1 Introduction 

Teleparallel equivalent of General Relativity (TEGR) [1, 2] is an equivalent formulation of 
gravity, but, instead of using curvature invariants defined by the Levi-Civita connection, 
the Weitzenbdck connection is adopted. Therefore TEGR exhibits no curvature but only 
torsion. The dynamical objects in such a framework are the four linearly independent 
vierbeins and the advantage of this framework is that the torsion tensor is formed solely by 
products of first derivatives of these vierbeins. In such a formulation, as described in [2], 
the Lagrangian density, T, can then be constructed from this torsion tensor assuming the 
invariance under general coordinate transformations, global Lorentz transformations, and 
the parity operation, along with demanding the Lagrangian density to be second order in 
the torsion tensor. In [3-6] an extension of the above idea was constructed, making the 
Lagrangian density a function of T, similar to the f(R) extension of the Hilbert-Einstein 
action. f(T) gravity has gained a significant attention in the literature, and proves to 
exhibit interesting cosmological implications [3-51]. 

Such an approach can be framed within the class of new gravity theories aimed to 
extend General Relativity in order to solve its shortcomings at Infra- Red and Ultra- Violet 
scales [52]. Clearly, in extending the geometry sector, one of the goals is to solve the 
puzzle of dark energy and dark matter that, up to now, seems to have no counterpart at 
fundamental level. In other words, both f(T) gravity and f(R) gravity could be reliable 
approaches to address the problems of missing matter and accelerated expansion without 
asking for new material ingredients that have not been detected yet by the experiments 



[53]. 
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In this work we investigate D-dimensional f(T) gravity, considering additionally the 
electromagnetic sector. Exact black-hole solutions with flat transverse section (Banados- 
Teitelboim-Zanelli (BTZ)-like solutions [54]) can be derived for a given range of parameter 
space. Then we analyze the singularities of these solutions based on the torsion scalar and 
the curvature scalar, pointing out differences with respect to f(R) gravity. It is important 
to stress that searching for exact solutions is a fundamental step to set a new field theory. 
Exact solutions allow a full control of the systems and can contribute to the well- formulation 
and well-position of the Cauchy problem (for a discussion on this point see [55]). 

The paper is organized as follows. In Sec. 2, we present a brief review of TEGR in 
four dimensions, as well as of its f(T) extension. In Sec. 3 the D-dimensional teleparallel 
gravity is formulated and the analysis is extended to D-dimensional Maxwell-/ (T) gravity. 
In Sec. 4 we derive exact charged static solutions and Sec. 5 is devoted to the investigation 
of singularities and horizons. Finally, in Sec. 6 we discuss some physical implications of 
the results. 

2 Teleparallel Equivalent to General Relativity and its f(T) extension 

In this section we briefly review Teleparallel Equivalent to General Relativity (TEGR) in 
four dimensions and its f(T) extension. Throughout the manuscript we use the following 
notation: Greek indices [i,v,... run over all space-time coordinates 0, 1, 2, 3; lower case 
Latin indices (from the middle of the alphabet) ... run over spatial coordinates 1, 2, 
3; capital Latin indices A,B,... run over the tangent space-time 0, 1, 2, 3, and lower case 
Latin indices (from the beginning of the alphabet) a, run over the tangent space spatial 
coordinates 1, 2, 3. 

As we mentioned above, the dynamical variable of teleparallel gravity is the vierbein 
field e^(x /i ), which forms an orthonormal basis for the tangent space at each point x^ 
of the manifold, that is • e# = T]ab, with t]ab — diag(l, — 1, — 1, — 1). Moreover, the 
vector &a can be expressed in terms of its components in a coordinate basis, namely 
e A = e A^M- I n sucn a formulation the metric tensor is acquired from the dual vierbein as 

9^u{x) = t]ab e p (x) ef (x) . (2.1) 

Although in General Relativity one uses the torsionless Levi-Civita connection, in the 
present construction one uses the curvatureless Weitzenbock connection [56] , whose torsion 
tensor reads 

W A W A 

T\ v = - IV = e\ (d^ - d v e£) . (2.2) 

Moreover, the contorsion tensor, which gives the difference between Weitzenbock and Levi- 
Civita connections, is given by K^ v p = — \ (T^p — T ViL p — T p while it proves convenient 
to define Sp^ = \ (K^ v p + 5 P T au a — 5 p T a „) . For a detailed exposition of torsion tensor 
properties see [57]. 

In conclusion, in the present formulation the torsion tensor includes all the infor- 
mation concerning the gravitational field. Using the above definitions one can construct the 
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simplest form of the "teleparallel" Lagrangian, which is the torsion scalar, that is [58, 59] 
C = T = -JP^T^ + \t^T vw - T pf fT^ u . (2.3) 
Thus, the simplest action of teleparallel gravity reads: 

S =±- f d 4 xe(T + C m ) , (2.4) 

where k = 8irG, e = det(e^) = \J—g and C m accounts for the matter Lagrangian. It is 
worth noticing here that the Ricci scalar R and the torsion scalar T differ only by a total 
derivative of the torsion tensor, namely [60]: 

R=-T-2V»{V I J). (2.5) 

Varying the action (2.4) with respect to the vierbein we obtain the field equations 

1 em 

e-%(ee? A S p n - e A T% x S p ^ - -e A T = ^Ge p A T / , (2.6) 

em 

where the tensor T p v on the right-hand side is the usual energy-momentum tensor of 
matter fields. These equations coincide with those of General Relativity for every geom- 
etry choice, and this is the why the theory is named "Teleparallel Equivalent to General 
Relativity" . 

One can generalize the above formulation considering arbitrary functions of the torsion 
scalar f(T) in the gravitational action [3-6], although the Lorentz invariance of the linear 
theory seems to be spoiled [61, 62]. Thus, the action becomes 

s = hj d * xe [T + /(T) + Cm] ■ (2,7) 

Notice the difference in the various conventions in 4D-/(T) literature, since some authors 
replace T by f(T), while the majority replace T by T + f(T). In this work we follow the 
second convention, that is teleparallel gravity is acquired by setting f(T) = 0. Finally, 
variation of the action (2.7) with respect to the vierbein gives the field equations 

/ (if \ fffl em 

e-%(ee A S p n \l + -e A T^ x S p ^+e A Srd,(T)^--e A [T+f(T)} = 4irGe A T /. 

(2.8) 

3 D-dimensional Teleparallel Gravity and its Maxwell-/ (T) extension 

In this section we present teleparallel gravity in D-dimensions and its Maxwell-/ (T) ex- 
tension and we explore its properties. It proves more convenient to use differential forms, 
where the torsion 2-form T a is simply T a = de a . 

We start with the gravitational teleparallel action with the most general quadratic 
form in the torsion tensor. Under the assumption of zero spin-connection it is given by 
[63, 64] 



— / (p C + p\Ci + p 2 C 2 + P3-C3 + , (3.1) 

\k J 
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where pi are dimensionless parameters and 

£ = ^e a A *e a , L\ = de a A *de a , £ 2 = (de a A *e a ) A *(cfe& A e b ) , 

£ 3 = (de a A e b ) A *(de a A e b ) , £4 = (de a A *e b ) A *(de 6 A e a ) , (3.2) 

with * standing for the Hodge dual operator and A for the usual wedge product. The 
coupling constant po = — |A accounts for the cosmological constant term, and furthermore, 
since £3 can be completely expressed in terms of C\, in the following we set /O3 = 
[63]. Lastly, we mention that in the above expression k is the D-dimensional gravitational 
constant, while the vierbeins and the metric are now D-dimensional. Therefore, in the 
following, all the conventions adopted in Sec. 2 extend in D dimensions. 
Action (3.1) can be written more conveniently as 

S= -J- /"(T-2AW1 = -!- [ d D xe(T -2A), (3.3) 
2k, J 2k J 

where *1 = (— l) D ~ 1 e° A e 1 A e 2 A e D ~ l , and the torsion scalar T is given by 



T = (-1) * pi(de a A *de a ) + p 2 (de a A e a ) A *(de b A e b ) + p 4 (de a A e 6 ) A *(de b A e a ) 

(3.4) 

Expanding this expression in its components we acquire 

T = \ (pi + p 2 + Pi) T abc T abc + P2 T abc T bca - PA T^T b c , (3.5) 

thus we straightforwardly see that for p\ = 0, p 2 = — \ and p\ = 1 it coincides with (2.3) 
in D dimensions, namely 

T = \T abc T abc - X -T abc T bca - T a ac T b bc . (3.6) 

Now, we will extend the above discussion considering arbitrary functions of the torsion 
scalar f(T) in the D-dimensional gravitational action. Thus, we consider an action of the 
form 

s = hj dZ M T + /( T )- 2A ] . ( 3 - 7 ) 

with the torsion scalar T given by (3.5), that is we keep the general coefficients pi. In 
differential forms the above action can be written as 

S = Yk\ {[/( T ) + T - 2A ]*!} > ( 3 - 8 ) 

where now T is given by (3.4). Finally, note that teleparallel D-dimensional gravity dis- 
cussed above is obtained by setting f(T) = 0. 

Lastly, we extend the discussion incorporating additionally the electromagnetic sector. 
In particular, we extend the total action to 

S= TkI {[/( T ) + T - 2A ]*!} + / £f, (3.9) 
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where 



--F A* F 
2 



(3.10) 



is the Maxwell Lagrangian, while F = dA, with A = A^dx^, is the electromagnetic poten- 
tial 1-form. The action variation leads to the following field equations: 

5C = ^A{(l + jpi 2d-k de a + i a (de b A *de b ) - 2i a (de b ) A *de b 

+P2 |— 2e Ad* (de b A e b ) + 2de a A *(de b A e b ) + z a de c A e c A *(cfe b A e&) 

-2i a (de b ) Ae b A *(de c A e c )| 
+/?4 |— 2e b A d* (e a A de b ) + 2de b A *(e a A de b ) 

+i a e c A de b A *(de c A e b ) - 2i a (de b ) A e c A *(de c A e fe )| j 



d 2 f 



pi * de a + p 2 e a A *(de b A e b ) + p±e b A *(de b A e a ) 



+ 



A *e a - 2A * e a - - [F A i a (*F) - i a (F) A *F] 



+5 A (d*F) = , 



(3.11) 



where i a is the interior product. Although one could investigate solution subclasses with 
general coupling parameters pi, in the following, for the sake of simplicity, we restrict to 
the standard case p\ = 0, p 2 = —1/2 and p± = 1 of (3.6). 



4 Exact charged solutions 

Let us now investigate the charged solutions of the theory. In order to extract the static 
solutions we consider the metric form 

i=D-2 

ds 2 = F {rf dt 2 —dr 2 ~r 2 dx f > 

G ( r ) i=i 

which arises from the vierbein diagonal ansatz 

e° = F (r) dt , e 1 = -prT^ dr > ^ = rdxi > ^ = rdx2 ' ( 4 - 2 ) 

G (r) 

Let us make an important comment here concerning the vierbein choice that corresponds 
to the metric (4.1). In the case of linear-in-T gravity the above simple, diagonal relation 
between the metric (4.1) and the vierbeins (4.2) is always allowed. On the contrary, in 
the extension of f(T) gravity, in general, one could have a more complicated relation 
connecting the vierbein with the metric, with the vierbein being non-diagonal even for a 
diagonal metric [65-68]. However, in the cosmological investigations of f(T) gravity [3— 
51], as well as in its black-hole solutions [69-73], the authors still use the diagonal relation 
between the vierbeins and the metric, as a first approach to reveal the features of the 
theory. Thus, in the present investigation we also impose the diagonal relation between 
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the vierbeins and the metric, as a first approach on the subject and in order to reveal the 
main features of the solution structure. However, we are aware that a detailed study of the 
general vierbein choice (and its relation to extra degrees of freedom) is a necessary step for 
a deeper understanding of /(T)-gravity foundations. 

Concerning the electric sector of the electromagnetic 2-form we assume 



i=D-2 



dA 



E r {r)e 1 Ae° + ^ Ei(r)e°Ae 



t+l 



(4.3) 



i=l 



where E r is the radial electric field, neglecting for the moment the magnetic part. Thus, 
inserting the above ansatzes in the field equations (3.11), we finally obtain 



1 + 



dT 



dT 



i=D-2 



. 



G (r) G' (r) + F' (r)G(rY 



rF (r) 



~ 2 



i=l 



(4.4) 



^'(^^-f^E^O, (4.5) 



1 + 



df_ 
dT 



F" (r) G {rf F' (r) G' (r) G (r) F' (r) G (r)' 



dT 2 



T' (r 



F(r) 



1 + — 



F' (r)G(rf 



(79-3) 



F(r) 
G'(r)G(r) 



rF (r) 



(S-3) 



F(r) 



+ (D - 3) 



where 



T(r) 



2{D 



= i = 1,. ..,£>- 2 , 
i,j = l,...,D-2 (i^j) 



F' (r) G (r) 



+ (D — 2) (D — 3) 



G(r) 



(4.6) 

(4.7) 
(4.8) 



(4.9) 
j equation 



rF (r) 

The remaining field equations are equivalent to equation (4.6), that is the a 
is similar to (4.6), but with —\E\ replaced by —\E 2 -_ X . 

A first observation is that from (4.4) we deduce that T has, in general, an r-dependence, 
which disappears for a zero electric charge. Such a behavior reveals the new features 
that are brought in by the richer structure of the addition of the electromagnetic sector. 
Moreover, form (4.7) and (4.8), we deduce that we cannot have simultaneously two non-zero 
electric field components. This result is similar to the known no-go theorem of 3D GR-like 
gravity [74, 75], which states that configurations with two non- vanishing components of 
the Maxwell field are dynamically not allowed. However, it is not valid anymore if we add 
the magnetic sector, as we will see in subsection 4.3 (it holds only for D=3). Therefore, 
in the following we investigate the cases of radial electric field, of non-radial electric field, 
and of magnetic and radial electric field, separately. 
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4.1 Radial electric field 

We first consider the case where there exists only radial electric field. Thus, the Maxwell 
equations give 

E r = -§^, (4.10) 



where Q is an integration constant which, as usual, coincides with the electric charge of 
the black hole. Now, integrating equation (4.5) we find the very simple and helpful result 



F(r) = G(r)[l + 



dT 



(4.11) 



Using equations (4.11) and (4.9) we obtain 



dG {rf 
dr 

whose solution is 
G(r) 2 = 



+ 



dr 



In 1 + 



dT 



+ 



(D-3) 



G(r) 



rT (r) 



1 + J£ 

1 + dT 



„D-3 



(D 



1 + 



4f 
dT 



„D-2 



T (r) dr + Const 



(4.12) 



(4.13) 



and using equation (4.11) we get 



F(r) 



,D-Z 



(D 



1 + ) r D ~ 2 T 



(4.14) 



where Const is an integration constant related to the mass of the spherical object. 

In order to proceed, and similar to [73], we will consider Ultraviolet (UV) corrections 
of f(T) gravity. In particular, we examine the modifications on the solutions caused by UV 
modifications of D-dimensional gravity and we consider a representative ansatz of the form 
f(T) = aT 2 . This is the first order correction in every realistic f(T) gravity, in which we 
expect f(T) <^ T [14, 40], since T (like R) is small in K 2 -units. Thus, for q/0, equation 
(4.4) leads to 

-1 ± v 7 ! - 24aA - 6aQ 2 r 4 - 2D 



T(r) 



6a- 



(4.15) 



with the upper and lower signs corresponding to the positive and negative branch solutions 
respectively (note that if a = then (4.4) becomes linear having only one solution, which is 
given by the a — > limit of the positive branch of (4.15), namely T(r) = —Q 2 r A ~ 2D /2 — 2A, 
in which case teleparallel gravity is restored). Thus, 



l + ^ = -±-v/r- 24qA - 6aQ 2 r 4 " 2C , 
dT 3 3 V ^ 



(4.16) 
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and therefore performing the integration that appears in (4.13) and (4.14), we obtain 



1 + 



dT 



r D ~ 2 T(r)dr 



18aQV 



yJr iD (1 - 24aA - 6aQ 2 r 4 ~ 2D ) 
54a 



(1 + 72aA) r D ~ l 
3-D D - 1 

6aQ 2 r 3 ~ 3D (-1 + 2AaA) 



2D - 5 



D - 1 



(D - 2f (-1 + 24aA) Q 2 r 



1 + 



6aQ 2 r 4 - 2P ^ / D-3 1 3D-7 . 6aQ 2 r 4 - 2D 
2ri \ 2(D-2) ' 2' 2(D-2) ' l-24aA 



-1+24«A 



,(4.17) 



3 (D - 3) (2D - 5) (D - 1) y/r 4D (1 - 24aA - 6aQ 2 r 4 ~ 2D ) 



where 2-^1 ( a ) &j c; x) is the hyper geometric function. We mention that the last argument of 
this function, namely (6aQ 2 r 4 ~ 2D ) / (1 — 24aA), must be negative, while from (4.15) it is 
required that 1 — 24aA — 6aQ 2 r 4 ~ 2D must be positive, therefore we deduce that a should 
be negative. 

In summary, inserting the integral (4.17), along with (4.16), in (4.13) and (4.14), we 
find that the black-hole solution is: 



G{rf 



54a 



3 - 1 - 3 

18aQ 2 r 3 ~ D 



\(D- 



l ± iy 7 ! " 24aA - 6aQ 2 r 4 " 2B ) r D-3 1 C° " 2 ) 
(1 + 72aA) r 13 - 1 " 



D 



D - 1 



± 



y/r 4D (1 - 24aA - 6aQ 2 r 4 ~ 2D ) r6aQ 2 r 3 ~ 3D (-1 + 24aA) r' 1 - 



54a 

(D-2) z (-l + 24aA) Q 2 r 



2D 



^2,,3+D. 1 1 + 6aQ 2 r 4 ~ 2 



f+24aA 21 I 2(D-2)> 2' 2(D-2) ' l-24aA 



D - 1 

D-3 1 3D-7 



3 (I? - 3) (2D - 5) (D - 1) vV 4D (1 - 24aA - 6aQ 2 r 4 ~ 2D ) 



+Const 



(4.18) 



and 



F(r) 



— i — 

D -a I CD — 



r D-3 I (D - 2) 

lSa^r 3 -- (1 + 72aA) r 13 " 1 



54a 



3-D 



D - 1 



± 



^ r 4D (1 _ 24 a A - 6aQ 2 r 4 ~ 2D ) r6aQ 2 r 3 - 3D (-1 + 24aA) r" 1 - 



54a 

(D-2) z (-l + 24aA) Q 2 r 



D 



2D - 5 



D - 1 



2, 1 n^n^+O./i.fegV^ f 

x 1 1 I O/l^,, A Z 1 1 



D-3 1 3D-7 . 6aQ 2 r 4 - 2P 
l+24aA 2^1 I 2(Z?-2)' 2' 2(D-2) ' l-24aA 



3 (D-3) (2D - 5) (D - 1) yV 40 (1 - 24aA - 6aQ 2 r 4 ~ 2D ) > 
+Consty (4.19) 

The special point in the parameter space A = l/(24a) needs to be analyzed separately, 
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since in this point we obtain the solution 



1 + f;) 



-1-5D 



27a 



2r 6D 9aQ 2 r 4 + 4D 3^6 (-aQ 2 r 4+2D ) 3/2 



D - 1 



D 



2D 



(4.20) 



and thus 



G(rY 



| ± \^/l-24aA-6aQ 2 r i ~ 2D ) 2 r D ~ 3 1 (jD " 2) 



-1-5D 



2r 6D 9aQ 2 r 4+4D 3^6 (-aQ 2 r 



2^4+2D\3/2 



27a ID — 1 



D 



2D 



}+ Const}, (4.21) 



and 



F(rf 



1 



1 



D ~ 3 I (D — 2) 



-1-5D 



27a 



2r 6D 9aQ 2 r 4+4D 3^ (-aQ 2 r 4 + 2Z) ) 3/2 



D — 1 



2D 



}+ Const}. (4.22) 



Finally, the case -D = 3 has to be analyzed separately. Taking properly the limit 
D = 3 of the above expressions we obtain the solutions extracted in [73] for 3D Maxwell- 
f(T) gravity. Lastly, one can straightforwardly check that in the limit a — > (of the 
positive branch since in this case the negative branch disappears) one re-obtains the usual 
charged General Relativity solutions. 

4.2 Zero radial field 

Let us for the moment assume that we have zero radial field. In this case equation (4.8) 
implies that we can have at most one non-zero component of the electric field along the 
non-radial (transversal) directions. However, as we mentioned below equation (4.9), for 
D > 3 the remaining field equations are similar to equation (4.6) but with — \E\ replaced 
by —^Ej_ 1: therefore subtracting these equations we acquire the conditions Ef = E 2 -, with 
% and j running from 1 to D — 2. These conditions, along with equation (4.8), yield E{ = 
(i = 1, ...,D — 2) for D > 3, that is the electric field is completely zero. The only cases 
where zero radial electric field does not lead to a disappearance of the total electric field 
is for D = 3 (where a non-zero azimuthal electric field is possible) which was analyzed 
in detail in [73], or if we consider simultaneously non-zero non-radial electric field with 
magnetic field, case which lies beyond the scope of the present investigation. 

4.3 Magnetic field and radial electric field 

For completeness we also examine the case where magnetic field is present. While in D = 3 
we deduce that electric field must be absent [73], for D > 3 one can simultaneously have 
non-zero magnetic and electric fields. As an explicit example we consider an electromag- 
netic strength 2-form in four dimensions given by 



E r (r) e 1 A e° + B 23 (r) e 2 A e 3 



(4.23) 
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that is we consider a radial electric field E r and a magnetic field B23 both depending on the 
radial coordinate r only. Prom the Maxwell equations in four dimensions for the electric 
field we immediately obtain 

E r (r) = Q , (4.24) 

while incorporating the equations of motion analogous to (4.4)-(4.8) we can see that a 
solution is obtained by 

£ 23 (r) = ^, (4.25) 

leading to the metric coefficients (4.19) and (4.20) with Q 2 + P 2 in place of Q 2 (and for 
D = 4). 

5 Singularities and horizons 

Let us now investigate the singularities and the horizons of the above solutions. The 
first step is to find at which r do the functions G (r) 2 and F (r) 2 become zero or infinity. 
However, since these singularities may correspond to coordinate singularities, the usual 
procedure is to investigate various invariants, since if these invariants diverge at one point 
they will do that independently of the specific coordinate basis, and thus the corresponding 
point is a physical singularity (note that the opposite is not true, that is the finiteness of an 
invariant is not a proof that there is not a physical singularity there) . In standard black- hole 
literature of curvature- formulated gravity (either General Relativity or its modifications), 
one usually studies the Ricci scalar, the Kretschmann scalar, or other invariants constructed 
by the Riemann tensor and its contractions. 

In teleparallel description of gravity, one has, in principle, two approaches of finding 
invariants. The first is to use the solution for the vierbein and the Weitzenbock connection 
in order to calculate torsion invariants such as the torsion scalar T. The second is to use the 
solution for the corresponding metric in order to construct the Levi-Civita connection, and 
then use it to calculate curvature invariants such are the Ricci and Kretschmann scalars 
(a calculation of curvature scalars using straightaway the Weitzenbock connection leads to 
zero by construction). The comparison of both approaches is a main subject of interest of 
the present work, capable of pointing out differences between curvature and torsion gravity. 

In particular, we are going to investigate whether one can formulate everything in 
terms of vierbens, Weitzenbock's connection and torsion invariants, as one can do with the 
metric, the Levi-Civita connection and curvature invariants. Perhaps one could say that 
the use of curvature invariants, instead of torsion ones, is better justified by the fact that in 
a realistic theory matter is coupled to the gravitational sector through the metric and not 
through the vierbeins (with the interesting exception of fermionic matter), and particles 
follow geodesies defined by the Levi-Civita connection. On the other hand, one could say 
that the two approaches are equivalent only with a suitable, non-diagonal, relation between 
the vierbeins and the metric. In any case, while at the classical level the above approaches 
could look equivalent or alternative, for the quantization procedure it would be crucial 
to determine whether the metric or the vierbein is the fundamental field. Therefore, the 
following analysis can enlighten this subject. 
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The torsion invariant T, that is the torsion scalar, that arises from the vierbein solution 
(4.2) with the use of Weitzenbock's connection is (4.9), which in the examined case becomes 
just (4.15). The curvature invariants that arise from the metric solution (4.1) through the 
calculation of the Levi-Civita connection are 



= G{rfF"{r) _ G (r) G> (r) F> (r) _ 2(D _ 2) G(r) 2 F'(r) 
F(r) F(r) K ' rF (r) 

_ 2 (D - 2) G(r)G/(r) - (D - 2) (D - 3) , (5.1) 



R ^ RlLV = G } r l \ V F " M G M + rF ' (0 G ' W + (d-2)g (r) f' (r)] 2 

F(r) r 2 

+ _GyO_ r F // ^ G ^ + rF , ^ qi ^ + (£, _ 2 ) G' (r) F (r)l 2 
F (r) r 2 

+ (D - 2) [rG (r) F' (r) + rG' (r) F (r) + (D - 3) G (r) F (r)] 2 ,(5.2) 

RrpaBT"" = 4^£ [F" (r) G (r) + F' (r) G' (r)] 2 + 4 (D - 2) G ®- 1 ' ^ 

+4 (£> - 2) G/ (r) 'f (r)2 + 2 (J - 2) (D - 3) ^ , (5.3) 

being respectively the Ricci scalar, the Ricci tensor square and the Kretschmann scalar. 
Note that, using (4.9), the Ricci scalar is given by 

R = _ T - 2 G{r)2F " {r) - 2 G M&MI1M - 2 (D-2) G ^ G '^ , 
F (r) F (r) r 

which is just relation (2.5) calculated for the vierbeins (4.2). 

Observing the form of the torsion scalar T in (4.15) we deduce that in the charged 
case it diverges only at r = 0. This can be alternatively verified examining the form (4.9) 
along with the expressions (4. 18), (4. 19) (or (4. 21), (4. 22) for the special solution). 

Observing the forms of Ricci and Kretschmann scalars in (5.1), (5.3) we deduce that 
in the charged case the possible divergence points are at r = 0, at the points where 
G(r) 2 —7- oo, or at the roots of F(r). From the solutions for G(r) and F(r) of (4. 18), (4. 19) 
we straightforwardly obtain that r = indeed leads to divergent Ricci and Kretschmann 
scalars. From the form of G(r) in (4.18) along with (4.16) we observe that G(r) 2 — > oo at 
1 + |L = 0, that is at 

1 + 8aA\ 4=2£> 



a relation which holds only for the negative branch, since for the positive branch 1 + ^ has 
no roots (additionally since a < the above solution is real only for A < — l/(8a)). Indeed 
one can straightforwardly see that the Ricci and Kretschmann scalars do diverge at r = r s . 
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Finally, concerning the roots of F(r), due to (4.11), namely F (r) = G (r) 



1 + Trr J i the y 



are just the roots of G(r), since F(r) remains finite and non-zero at 1 + -Jp = since at 
this point G(r) 2 — > oo. Taking the corresponding limits and using (4. 18), (4. 19), we can 
see that the roots of G(r) leads always to finite Ricci and Kretschmann scalars. All the 
above hold also for the special solution (4. 21), (4. 22). In summary, in the charged case the 
Ricci and Kretschmann scalars diverge at r = and at r = r s given by (5.4). Lastly, we 
mention that in the uncharged case the Ricci scalar vanishes, however the Kretschmann 
scalar behaves as in the charged case. 

From the above analysis we are led to the very interesting result that in some cases 
the singularities obtained by the torsion scalar analysis are less than those obtained by 
the curvature scalar analysis. In particular, this happens for the negative branch of the 
solutions, for Q / and for A < — l/(8a), in which case the curvature invariants possess an 
additional physical singularity at r = r s given by (5.4). We stress that when f(T) = 0, that 
is in the case of usual teleparallel gravity, the negative branch disappears as we mentioned 
above, thus the singularity analyses of the two approaches coincide. Additionally, in the 
uncharged case, that is when Q = 0, the extra singularity at r s disappears too, and the 
singularity analyses of the two approaches coincide too. In conclusion, we deduce that 
the above difference in the physical singularities of the torsion and curvature analysis, is 
a result of both the non-linear f(T) structure and of the non-zero electric charge, which 
reveals the novel features that are brought in in the theory in this case. 

Let us discuss on the horizons of the above solutions. Although we showed that the 
roots of G(r) at r > (if they exist) do not correspond to physical singularities, obviously 
they correspond to horizons, since they appear in the denominator in the metric (4.1) 
and in the vierbeins (4.2). In order to show that the roots of G(r), say at r = r#, 
are just coordinate singularities, we consider the Painleve-Gullstrand coordinates [76-79] 
through the transformation dt = dr+g (r) dr, with g (r) a function of the radial coordinate. 
Therefore, the metric (4.1) becomes 



ds z 



F (r) 2 dr 2 + 2g (r) F (rf drdr 



G(r) 



Firfgirf 



dr' 



=D-2 



Choosing g (r) 



F{rY 

h(rf = 



Girf 



and defining 



F{rf /G{rf 



2 + v 7 ! - 24aA - 6aQ 2 r 4 " 2D /3 



we can bring it in a flat Euclidean form 



(5.6) 



i=D-2 



ds 2 = F (rf dr 2 + 2h (r) \ll-G (rfdrdr - dr 2 - r 2 



dx 2 , 



i=i 



(5.7) 



which is regular at r = rn- Therefore, r = r#, if they exist, are just coordinate singularities, 
that is horizons. 

From the above analysis it is implied that the black-hole solutions of the charged f{T) 
gravity may possess a horizon at r# that shields the physical singularities. However, firstly 
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it is not guaranteed that ru exists, since there could be parameter choices for which G(r) 
has no roots, that is the physical singularity at r = becomes naked. Secondly, even 
if ru exists it is not guaranteed that it will shield the second physical singularity of the 
charged negative branch at r = r s given by (5.4), since this will depend on the specific 
parameter choice. In particular, we can see that if F(r s ) 2 < then r# exists and shields 
the singularity at r s , that is ru > r s , otherwise r s is a naked singularity. This is not the 
case for f(T) — > or Q — > 0, in which, as we mentioned, r s disappears. Therefore, we 
conclude that the cosmic censorship theorem, namely that there are always horizons that 
shield the physical singularities, does not always hold for f(T) gravity, a result that was 
already found in the 3D case too [73]. 

Before proceeding to the numerical elaboration of the obtained solutions, we make 
the following comment. In curvature gravity there can be cases where the Ricci scalar is 
finite at one point although there is a physical singularity there, which is revealed through 
the use of the Kretschmann scalar, and that is why people usually examine both scalars 
simultaneously. Thus, one could ask whether one should use additional torsion scalars too, 
defined as various contractions of the torsion tensor. In particular, according to (2.3), the 
torsion scalar T contains three separate scalars, corresponding to different contractions of 
the torsion tensor, namely I\ = T^PT^p, I2 = T pfJ,u T Ufip , I3 = T PfjL P T U ' J ' U , and thus one 
could additionally examine their behavior in order to reveal the singularities. However, 
it is well known that these separate combinations are not invariant under local Lorentz 
transformations, and that was the reason that the teleparallel Lagrangian (torsion scalar) 
T was defined as their specific combination which becomes Lorentz invariant [2]. Therefore, 
one cannot use other torsion scalars apart from T in order to investigate the singularities 
(the explicit calculation of Ij's for the obtained vierbein solutions shows that they acquire 
different values in different coordinates, and thus they are not invariants, however their 
specific combination in T does acquire the same value independently of the coordinate 
basis and thus it is a well-defined invariant). 

In order to provide a more transparent picture of the above singularity and horizon be- 
havior, we proceed to the numerical elaboration of specific examples. Since the asymptotic 
behavior of G(r) is given by 

G(r) 2 = -A e// r 2 + ... , (5.8) 
where . . . correspond to sub-leading terms and 

A eff = ; [l + 72aA =f (1 — 24aA) 3 / 2 

t} 54a(L> - 1)(D - 2)(2 ± y/l - 24aA) L ^ V ' 

we can distinguish three subclasses, namely the asymptotically AdS one (A e jf < 0), the 
asymptotically dS one (A e jy > 0) and the limiting A. e ff — one. Note that the second and 
third solution subclasses can never be obtained by the negative solution branch. Without 
loss of generality we consider D = 4, a = —1 and Const = —1, while we suitably choose A 
in order to lie in the above three subclasses, which we investigate separately. 



(5.9) 
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• Case A e fj < 

Let us consider A = —1/25, which satisfies the condition A e ff < 0. In Fig. 1 from up 
to down we depict G(r) 2 , the Ricci scalar, the Kretschmann scalar and the torsion 
scalar T{r) as a function of r, for the positive and negative branches of black-hole 
solutions. The left graphs correspond to charged solutions (Q = 1), while the right 
graphs correspond to uncharged solutions (Q = 0). 

The positive branch in the charged case exhibits only one physical singularity at 
r = 0, in which the torsion scalar and the Ricci and Kretschmann scalars diverge, 
however it is shielded by two horizons at rjj = T- and m = r +, since in this case 
G(r) 2 has two roots. In order to examine whether r+ is a Killing horizon we see that 
the timelike Killing vector of the metric is e^dn = dt 1 , with norm e^e^ = gu = F (r) 2 
which vanishes at r = r + . Inside the horizon the Killing vector field is spacelike, 
while outside it is timelike, and thus it corresponds to a null hypersurface. Finally, 
for the uncharged solutions we can see that the Ricci scalar and the torsion scalar are 
constants, while the Kretschmann scalar diverge at r = (similarly to usual General 
Relativity) . 

For the negative branch in the charged case, the torsion scalar possesses only one 
divergence, namely at r = 0, however the Ricci and Kretschmann scalars possess 
two divergence points, namely at r = and r = r s as described above. However, 
in this specific numerical example, both these physical singularities are shielded by 
the horizon at m > r s , in which all invariants remain regular. In particular, r# is 
a Killing horizon, corresponding to an event horizon since the Killing vector field is 
timelike outside the horizon and spacelike inside. Finally, for the uncharged solutions 
the Ricci and torsion scalars are constants, but the Kretschmann scalar diverges at 
r = 0. However, this physical singularity is shielded by the horizon at r = rjj in 
which G(r) 2 becomes zero. 

In summary, we indeed verify that in the charged case and for the negative branch 
the curvature invariants contain an extra divergence at r = r s , that does not appear 
in the torsion invariant, revealing the novel features of charged f(T) gravity. 

• Case K e ff > 

We consider A = 1/25, which satisfies the condition A e jy > for the positive branch 
(as we mentioned below (5.9) the negative branch cannot lead to A e jy > 0). In 
Fig. 2 from up to down we depict G(r) 2 , the Ricci scalar, the Kretschmann scalar 
and the torsion scalar T(r) as a function of r, for the positive branch of black-hole 
solutions. The left graphs correspond to charged solutions (Q = 1), while the right 
graphs correspond to uncharged solutions (Q = 0). 

The charged case possesses a physical singularity at r = 0, where the torsion scalar 
and the Ricci and Kretschmann scalars diverge, however it is shielded a horizon 

1 Note that none of the metric coefficients depends on time and thus the manifold has a timelike Killing 
vector dt, and similarly since none of the metric coefficients depends on x, there exist (D — 2) spacelike 
Killing vector fields d Xi . 
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Figure 1. The solutions for G(r) 2 of (4-13), for the Ricci scalar R(r) of (5.1 ), for the Kretschmann 
scalar R^u p(T R^ vpo "( r ) °f (5-3) an d f or the torsion scalar T(r) of (4-9), as a function of r, for the 
positive (thick solid curve) and negative (thin dashed curve) branch of the AdS solution subclass, 
for D = 4 ? a = — 1, Const = — 1 and A = —1/25. Left graphs correspond to charged solutions with 
Q = 1, while right graphs correspond to uncharged solutions with Q = 0. The thin 0-line is depicted 
for convenience. 
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Figure 2. The solutions for G(r) 2 of (4-13), for the Ricci scalar R(r) of (5.1), for the Kretschmann 
scalar R i i Vpa R ilv( ' a (r) of (5.3) and for the torsion scalar T(r) of (4-9), as a function of r, for the 
positive (thick solid curve) branch of the dS solution subclass, for D = 4, a = —1, Const = —1 and 
A = 1/25. Left graphs correspond to charged solutions with Q = 1, while right graphs correspond to 
uncharged solutions with Q = 0. The thin 0-line is depicted for convenience. 
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at r = rjj) where G(r) 2 becomes zero. In order to examine whether m is a Killing 
horizon, and similarly to the previous case, we observe that the timelike Killing vector 
of the metric is e^d^ = dt, with norm = gu = F (r) 2 which vanishes at r = r#. 
Since outside the horizon the Killing vector field is timelike, and inside it is spacelike, 
it is implied that it corresponds to a null hypersurface, that is a cosmological Killing 
horizon. 

In the uncharged case we observe that the Ricci scalar vanishes while the torsion 
scalar is constant, however the Kretschmann scalar diverges at r = 0. However, 
note that in this case G(r) 2 has no roots, that is there is not a horizon to shield the 
physical singularity at r = 0, which is therefore a naked one. 

• Case A e fj = 

We consider A = 0, which satisfies the condition h- e ff — for the positive branch 
(as we mentioned below (5.9) the negative branch cannot lead to A e ff = 0). In 
Fig. 3 from up to down we show G(r) 2 , the Ricci scalar, the Kretschmann scalar 
and the torsion scalar T(r) as a function of r, for the positive branch, with left 
graphs corresponding to charged solutions (Q = 1) and right graphs corresponding 
to uncharged solutions (Q = 0). 

The charged solutions possess a physical singularity at r = 0, where all invariants 
diverge, however it is shielded a horizon at r = ru, where G(r) 2 becomes zero. 
Examining the Killing vector, and similarly to the previous cases, we deduce that rn 
is a cosmological Killing horizon. In the uncharged case we observe that both the 
Ricci and torsion scalars vanish, however the Kretschmann scalar diverges at r = 0. 
Notice that in this case G(r) 2 has no roots, that is there is not a horizon to shield 
the physical singularity at r = 0, which is therefore a naked one. 

6 Concluding Remarks 

In this work we considered D-dimensional f(T) gravity including the Maxwell field. We 
extracted exact charged black- hole solutions depending on the functional form of f(T), on 
the electric charge and on the number the dimensionality D. Finally, we investigated the 
singularities and the horizons of the obtained solutions, following two different approaches. 
Firstly, by studying the torsion invariants constructed using the Weitzenbock's connection 
and the vierbein solutions, and secondly by studying the curvature invariants constructed 
using the Levi-Civita connection and the metric solutions. 

The main result is that in Maxwell-/ (T) gravity the curvature invariants possess more 
physical singularities than the torsion ones, in some particular solution subclasses. This 
difference disappears in the uncharged case, or in the case where f(T) gravity becomes the 
usual linear-in-T teleparallel gravity, thus it reveals the novel behavior that is introduced 
by the combined complication of the non-trivial f(T) structure with the electromagnetic 
sector. It seems that curvature and torsion invariants behave very differently depending 
on the presence of the Maxwell field. More generally, extending gravity in terms of f(T) or 
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Figure 3. The solutions for G(r) 2 of (4-13), for the Ricci scalar R(r) of (5.1), for the Kretschmann 
scalar R tlV p a W lvpc '(r) of (5.3) and for the torsion scalar T(r) of (4-9), as a function of r, for the 
positive (thick solid curve) branch, for D = 4, a = — 1, Const = — 1 and A = 0. Left graphs 
correspond to charged solutions with Q = 1, while right graphs correspond to uncharged solutions 
with Q = 0. The thin 0-line is depicted for convenience. 
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f(R) formulations could give very different results as soon as matter fields are taken into 
account. 

Finally, we have to note that, in the scenario we have considered, the physical singular- 
ities are not always shielded by horizons. Thus, the cosmic censorship does not always hold 
for D-dimensional Maxwell-/(T) gravity. From a cosmological point of view such a feature 
could be extremely relevant in order to investigate the early phases of cosmic evolution. 
On the other hand, considering astrophysical structures in strong field regimes, derived 
from torsion or curvature representation of gravity, could give rise to very deep differences 
in dynamics [22, 80]. 
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